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Abstract This series of three lectures covers (a) a basic introduction to 
symmetry breaking in general and chiral symmetry breaking in QCD, (b) an 
overview of the present status of lattice data and the knowledge that we have 
at finite temperature from chiral perturbation theory, (c) Results obtained 
from the Nambu-Jona-Lasinio model describing static mesonic properties are 
discussed as well as the bulk thermodynamic quantities. Divergences that are 
observed in the elastic quark-antiquark scattering cross-section, reminiscent 
of the phenomenon of critical opalescence in light scattering, is also discussed, 
(d) Finally, we deal with the realm of systems out of equilibrium, and examine 
the effects of a medium dependent condensate in a system of interacting 
quarks. 
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1 Introduction 



Chiral symmetry is the symmetry of Quantum Chromodynamics (QCD) that 
dictates the static properties of the low lying mesonic sector, in particular 
those pertaining to the pseudoscalar nonet (n,K,ri). This symmetry is re- 
sponsible for the fact that, in its broken phase, quarks acquire mass (and 
are termed "constituent" quarks, as they form parts of hadrons, while, in 
the restored phase, quarks have only their small or current mass values. It 
is believed that at finite temperature this symmetry is restored, a feature 
that is strongly motivated by numerical studies of QCD on the lattice. Con- 
comitantly with this picture, it is believed that another phase transition 
from a deconfined phase of matter (consisting of a hot fireball of quarks 
and gluons) to a confined phase can occur, in which only the final state of 
hadrons is observed. Given these two features, a large amount of scientific 
endeavor has been and will continue to be invested in the study of heavy 
ion collisions, in which high temperatures can be attained. In particular, the 
low-lying mesons are copiously produced, and since these provide the testing 
ground for chiral symmetry at T = 0, it is hoped that (with enough theo- 
retical and experimental study), a clear signal of this phase transition will 
emerge. To be quite precise, one requires unambiguous signals of both phase 
transitions, that of confinement /deconfinement, as well as chiral symmetry 
breaking/restoration. Thus far, however, there are no unambiguous signals 
known that are experimentally measurable for either of these transitions. In 
this paper, we shall confine ourselves to a discussion of chiral symmetry and 
its associated aspects, leaving the difficulties of confinement to a later stage. 

This series of three lectures is intended to introduce the concepts of chiral 
symmetry starting from basics. There is a short guide for the uninitiated into 
the ideas of what symmetry breaking is, and then an attempt to summarize 
the current status of what we know to be fact, taken on the theoretical 
level, at finite temperature. This involves examining firstly the lattice gauge 
simulations of QCD at finite temperature and then examining how far we 
can go with chiral perturbation theory From lattice gauge simulations, 

the existence of the chiral and deconfinement phase transitions is inferred. 
Critical exponents for the chiral transition have been obtained, but are as 
yet not conclusive. Temperature dependence of the mesonic screening masses 
have also been calculated, and the question of Ua(X) symmetry restoration 
addressed. Bulk thermodynamic properties have been studied over several 
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years, with larger and larger lattices, and this represents the state of the 
art of what we know today about these quantities in QCD. By contrast, 
while chiral perturbation theory gives a superb description of the low energy 
sector, and also gives the leading behavior expected of the order parameter as 
a function of temperature || || , it cannot per se be used to describe the phase 
transition region, which is non-analytic. The level of accuracy of CHPT at 
finite temperatures is illustrated in the calculation of the pion masses as a 
function of temperature in a recent publication || . 

Note that we restrict ourselves mainly to finite temperature and not to 
finite density. The first lattice simulations at finite density have already been 
performed . However, there are many technical difficulties that are not yet 
under control, and as such no results are completely reliable as yet. For 
this reason, we will also not attempt to make any model discussions at finite 
density at this stage, although there are of course several. 

In the second lecture, a simple chiral model, the Nambu-Jona-Lasinio 
(NJL) H H, [H| model is discussed, in which it becomes evident that fea- 



tures relating to static properties of the low-lying mesons are excellently 
reproduced. This includes charge radii, meson-meson scattering lengths, po- 
larizabilities, etc, and one can validate that the expected results of chiral 
perturbation theory are recovered, here with very few parameters. In ad- 
dition, the variation of the meson masses with temperature, although cal- 
culated in this model as pole masses, shows the same qualitative behavior 
as was observed by the lattice gauge groups. Given these successes with 
this model, one is encouraged to study the dependence of all static mesonic 
properties as a function of the temperature in order to investigate whether 
abrupt behavior occurs at the phase transition point. For two flavors of 
quarks, one finds that the pseudoscalar sector in particular is typified by an 
almost constant behavior in all static properties (such as the mean pion ra- 
dius (r 2 .) 1 / 2 , the polarizabilities a w and the scattering lengths a w ) for a wide 
range of the temperature shortly up until the point at which the chiral phase 
transition occurs, and then these quantities show a sharp divergence. This is 
true for the case in which the current quark mass of the up and down quarks 
mj} = wig = mo = 0, and a phase transition can occur. When itlq ^ 0, only a 
crossover can be observed in the order parameter. A new transition temper- 
ature T M = T Mott is defined as being the temperature at which the mesonic 
states become unbound, or resonances. It thus respresents a derealization 
of the mesons, rather than their deconfinement. The static properties of the 
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pionic sector then remain constant for most of the temperature range, and 
diverge at the Mott temperature. One thus still observes a dramatic struc- 
ture - either directly at the phase transition temperature itself in the case of 
m = or at T = T M for m 7^ 0. 

It is also extremely interesting to study dynamical quantities such as the 
elastic cross-section for qq — > qq scattering. This particular quantity displays 
a divergence at the critical or Mott temperature in a similar fashion as occurs 
in the phenomenon of critical opalescence that is observed in the scattering 
of light. However, although this feature and those observed for the static 
properties are exciting, their direct measurement is elusive if not downright 
impossible. 

The scalar mesonic sector within the NJL model is observed to display 
a completely different behavior. Here the mass drops relatively quickly with 
temperature. Nevertheless, experimentally, the scalar mesons constitute a 
multiplet that appears to have the symmetry badly broken, and the lowest 
meson of which (the a) has an extremely large width. Consequently only 
indirect information on this sector is useful. 

How then can one hope to observe the chiral phase transition? To at- 
tempt to answer this question, we recall that the chiral phase transition 
appears to be intimately linked with the confinement /deconfinement phase 
transition, i.e. they appear to take place at the same temperature fTI| . A 
heuristic understanding of this feature is quite satisfactory - it implies that 
at high temperatures, one should have chiral symmetry restored in a plasma 
phase, with free (current) quarks and gluons being the ingredients, while 
at T < T c , the confined phase contains only hadrons that are made up of 
constituent (massive) quarks. Experimental effort to detect the quark-gluon 
plasma phase is concentrated on contructing hot and dense matter via heavy- 
ion collisions such as Pb + Pb at increasingly high energies, and will form a 
main part of the program of the two accelerators RHIC at Brookhaven and 
the LHC (Geneva) that are currently under construction. Given the fact 
now that heavy-ion collisions take place over a small time scale, it is conceiv- 
able that the features of divergences occurring in both static and dynamical 
quantities might enter realistically into a non- equilibrium treatment of such 
collisions, which of course involves many particles, the lightest of which are 
the pions, and thus to measurable observables. 

For this reason, the final lecture is devoted to a discussion of non-equi- 
librium physics of an interacting fermionic Lagrangian, and which is then 



5 



applied to the Nambu-Jona-Lasinio model in the lowest possible terms in 
an appropriate double expansion in both h and the inverse number of colors 
N c [|T^, Using the simplest approximations that lead to a semi-classical 
result, one can recover a Boltzmann like equation for the quark distribution 
function. Here one sees that the problems are simply open ended. The issue 
of constructing interlinked equations dealing with several species of particle 
must be confronted and the issue of multiparticle production (hadronization) 
must be addressed, since the usual Boltzman collision scenario that incorpo- 
rates only binary collisions is inadequate for a relativistic description. 

Obviously it is an impossible task to discuss all aspects of chiral symmetry 
breaking and restoration within three lectures, and for this reason I have been 
highly selective in the material presented. There are many, many studies in 
the literature involving chiral symmetry, and I am in no way attempting 
in this paper to be comprehensive. The interested reader may also refer 
to the work of Refs.flT5[ for treatments of the linear sigma model at finite 



temperature, for example, and to the work of Ref. ||T5[ for discussions in the 
baryonic sector, in addition to the other general references that are given in 
the text. 

The structure of this manuscript reflects the three lectures directly: in 
Section 2, current factual information on the chiral transition, taken from 
lattice gauge simulations and chiral perturbation theory is presented. In 
Section 3, the Nambu-Jona-Lasinio (NJL) model is used to present the ram- 
ifications of symmetry breaking at the critical temperature. In Section 4, a 
non-equilibrium formulation of a theory of interacting fermions is described 
and the equations are investigated for the NJL model. In the concluding sec- 
tion, we discuss where this could possibly lead to observable consequences. 



2 Equilibrium thermodynamics. 

In this section, we attempt to present those aspects of chiral symmetry at fi- 
nite temperature that are regarded as being "exact" or factual, that is to say, 
they are derived from QCD itself, or from considerations thereof. We start by 
briefly introducing the reader to the general concept of symmetry breaking at 
T = 0. Following this, chiral symmetry breaking in the QCD Lagrangian is 
analysed. In the following subsection, the simulations of lattice gauge theory 
are discussed, dealing firstly with the temperature dependence of the order 
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parameter, the critical exponents obtained at the phase transition, meson 
screening masses and the question of whether Ua(1) symmetry is restored 
at high temperatures or not. Secondly, we indicate what is known from the 
lattice about bulk thermodynamic properties. The pressure density, energy 
density and entropy densities have been calculated on the lattice. These 
quantities give rather indications of the confinement /deconfinement transi- 
tion, and as we will show in Section 3, cannot be described well by a model 
that contains chiral symmetry alone, and which ignores the confinement as- 
pect. 

In the final subsection, we briefly introduce the concepts of chiral pertur- 
bation theory (CHPT) and we describe the state of the art results at finite 
temperature. As will be seen, these give an important functional depen- 
dence at low temperatures, but cannot be expected to cope with the phase 
transition region, which is non-analytic. 

2.1 Introduction to chiral symmetry at T = 0. 

The fact that a Hamiltonian, or equivalently a Lagrangian, is invariant under 
a symmetry transformation results in a degeneracy within the spectrum that 
is observed. Mathematically, one expresses the fact that a Hamilton function 
H is invariant under a specific symmetry via the statement 

UHtf = H (1) 

where U is an element of the group corresponding to this symmetry. Now if 
one considers the states \A) and \B) that are related by the transformation 

\B) = U\A), (2) 
it follows that \B) and \A) are degenerate, since 

E A = (A\H\A) = (B\H\B) = E B . (3) 

In order that this degeneracy manifest itself, however, it is necessary that 
the ground state of the system be invariant under such a transformation. 
Writing \A) and \B) in terms of creation operators, 

|A)=0 A |O) and |B)=0 B |O) (4) 
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with 

U<j) A tf = <Pb, (5) 
one sees that Eq.(0) holds only if 

\0) = U\Q), (6) 

i.e. the ground state is invariant under the symmetry group. Should this not 
be the case, one speaks of a spontaneously broken symmetry. 

Denoting U as U = exp(ie a Q a ) in terms of the (continuous) group pa- 
rameters e a and the generators of the symmetry 

Q a = J d 3 x,J«(x), (7) 

the statement Eq.(|6]) is seen to coincide with the equivalent form 

Q a \o) + o (8) 

although 

[Q a ,H]=0. (9) 

The direct consequence of this statement is that HQ a \0) = 0, or that there 
must exist a spectrum of massless particles with quantum numbers specified 
by the generators of the symmetry. This constitutes the Goldstone theo- 
rem. To be more precise, one can formulate this as follows: given that a 
Hamiltonian has continuous symmetries described by groups G\ requiring 
Nq 1 generators, while the ground state is invariant under groups G2 requir- 
ing Ng 2 < Nqx generators, the spontaneous breakdown of chiral symmetry 
leads to the existence of iV Gl — N G2 Goldstone bosons |16[ . 
Let us investigate now how this is applied to QCD. 



2.2 Chiral symmetry in QCD 

In this section, we analyse the symmetries of quantum chromodynamics, 
and compare this with the symmetry of the vacuum, determined purely by 
viewing the experimental spectrum. Start by examining the QCD Lagrangian 
itself, which can be written in a compact fashion as 

C QCD =$(iP- m )^ - \tr c G%G^ } (10) 
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where G^ u is the field strength tensor of the gluon field, 

G% = ~ dvG ^ — gfabcG^Gl, (11) 
is the co variant derivative, 

D^d^ + ig^G^x) (12) 

and f a bc are the structure constants of the symmetry group SU(3) RI7| . The 
quark field is a vector in flavor space, 

/ if) u {x) \ 
ip d (x 

i,= ^ (13) 



V ■ / 

and the (current) quark mass matrix is a diagonal matrix in flavor space, 

m Q = diag[m%,m%,mo, ...], (14) 

so that the second term in Eq.flTOD is 

irniQip = ^2m f tp f 'tp f . (15) 
/ 



If the quarks are massless, then the Lagrangian Eq. (|T0|) contains no term 
of the form Eq.flTB]) which can mix left and right handed components of the 
quark fields, that are defined as 



2(^7^, 



(16) 



i.e. these two fields are independent, and the Lagrangian remains invariant 
under transformations that individually transform these fields, 



^R,L - £W, U R , L E U{Nf), 



(17) 



and these are called chiral symmetries. However, a mass term of the form 
Eq.([H|) spoils this invariance since 



Terms ~ ipip = ipL^R + ^r^l 



(18) 
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mix left and right handed fields. Thus the term moipip constitutes an explicit 
symmetry breaking. 

The QCD Lagrangian Eq.fllUD is invariant under several transformations, 
such as 

V>->V' = e* aA > (19) 

etc. Accordingly, there are conserved Noether currents that correspond to 
these symmetries. They are 

Vg = t/^V A% = ip^-ip (20) 
V> = ^ = W75yV> (21) 

Among these is the current Aq, which corresponds to the transformation 
ip — > ip' = exp(z7 5 a:)?/>, where a is a continuous parameter. However, despite 
its appearance, this current is not conserved, 

3MJ = ^tr c G^. (22) 

This means that it does not reflect an underlying symmetry of the Lagrangian 
and its breaking was resolved by 't Hooft as being due to the presence of 



instantons 18 



One may thus identify the (continuous) symmetry groups of QCD as 
being generated by the charges of the remaining symmetries, and this is 

G x = SU L (N f ) ® SU R (N f ) ® U v (l). (23) 

On the other hand, by examining the particle spectrum that is observed 
experimentally, one finds that the symmetry of the vacuum is 

G 2 = SU V {N S )®U V {1). (24) 

Accordingly, there must be Nj — 1 massless Goldstone particles and these 
have the quantum numbers obtained from applying the axial charge operators 
to the vacuum, i.e. J p = O . In the case of two flavors, there are three such 
states, which are identified as corresponding to the charged and neutral pions. 
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For three flavors, one identifies the eight pseudoscalars as the pions, kaons 
and eta. One sees that the explicit symmetry breaking in this case is larger: 
wig ~ 150MeV in comparison with mjj ~ tuq ~ 5 MeV. 

The phase in which a system finds itself is usually characterized by an 
order parameter. This is an operator that transforms in a non-trivial fashion 
under the broken symmetry. Generally order parameters have the property 
of being zero in the symmetric or restored phase and non-zero in the sponta- 
neously broken phase, but this is not necessarily so. There are many possible 
ways of choosing an order parameter. The major criterion for doing so is 
that the order parameter should display that same invariances as the ground 
state. In the case of quantum chromodynamics, the ground state of QCD is 
invariant under Lorentz transformations and spatial reflections. The order 
parameter must thus be invariant under these same symmetries, and as such 
must be a scalar. The operator ipif) is the simplest choice. One thus makes 
the choice of (ipip), which is referred to as the quark condensate. 

2.3 Lattice gauge simulations 

Simulations of QCD on the lattice provide the most exact knowledge that 
we have of this theory that is derived from the QCD Lagrangian itself. The 
Lagrangian is discretized in space and time dimensions, and the variation 
with respect to the temperature of physical quantities is formally controlled 
by varying the size of the lattice in the temporal direction || f|], since 

T = l/N T a, (25) 

where a is the lattice size and iV r the temporal extent. In what follows, we 
simply list the major results that have been extracted via this methodology 
over the past few years. We show the temperature dependence of the chi- 
ral and deconfinement order parameters, discuss critical exponents, meson 
screening masses and {/a(1) resotration. Finally, we show plots of the bulk 
thermodynamic quantities. 

2.3.1 Order parameters 

The following recent results ^ [|] have emerged from the lattice gauge studies: 

• The pure gauge sector of QCD displays a well-established first order 
chiral transition at a rather high critical value of the temperature, T c = 
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270(5) MeV. The bulk properties for such a system are also well known 

El- ' 

• Full QCD including fermions displays a chiral phase transition at far 
lower critical temperature than that observed for pure gluonic systems. 
One finds T c ~ 150 MeV for two flavors of quark. 

• Studies of the Polyakov loop for quenched QCD places the critical tem- 
perature determined from the order parameter for deconfinement, To 
at about the same temperature at which the chiral transition T c occurs 

lj, i.e. 

Td — T c (26) 

This can be directly seen from Fig. 1, in which the order parameter for the 
chiral and deconfinement transitions are shown, together with their suscep- 
tibilities, as a function of j3 = Q/g 2 , over the transition region. Large (small) 
values of (3 represent the high (low) temperature regime. 

Based on these points, our physical (but heuristic!) understanding of the 
situation is that, at low energies, one has only hadronic states. These can be 
thought of as consisting of quarks carrying a dynamically generated quark 
mass m = m u = rrid for two flavors, and constructed into baryonic states 
or mesonic states according to the Goldstone theorem. At the temperature 
at which where chiral symmetry is restored, T c , and the constituent quarks 
take on their current mass value, deconfinement occurs simultaneously. The 
hadronic states dissolve, and one moves to a plasma containing only quark 
and gluonic degrees of freedom. 



2.3.2 Critical Exponents 

An obvious question that one may pose, when faced with a phase transition, 
is what are the critical exponents that govern the transition? Pisarski and 
Wilczek suggested that the dynamics of QCD is controlled by an effective 



scalar Lagrangian, constructed along the lines of the linear a model [|2T 



and, which for two flavors of quarks, has SU(2) <g> SU(2) = 0(4) symme- 
try. Now, according to arguments of universality [|22|, only the symmetry 
structure and dimensionality determine the values of the critical exponents, 
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m IT = . Of 




m /T = . Of 



Figure 1: Order parameters for the chiral and deconfinement transitions 
(lower and upper figures, respectively) are plotted as a function of the inverse 
QCD coupling (3 = 6/g 2 . High (low) values of (3 correspond to high (low) 
values of the temperature. The associated susceptibilities are also plotted in 



each case. Courtesy of [20 
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i.e. one expects that one should obtain the critical exponents of a 3D 0(4) 
symmetric spin model. The task of studying the critical exponents has been 
undertaken by a lattice group Q. Noting that the masses, which are respon- 
sible for explicit chiral symmetry breaking, play an analogous role to that of 
a magnetic field in the superconducting transition, these authors || £| have 
adopted the convention of defining a scaled quark mass as h — m q /T and the 
reduced temperature t = (T — T c )/T c . With this convention, the free energy 
density scales as 

f(t, h) = -^\nZ = b~ l f(b y % W*h), (27) 

introducing the thermal (yt) and magnetic (yu) critical exponents. Here b 
is an arbitrary scaling factor. There are various scaling relations that can 
be derived using Eq. (|27|) . In particular, one can show that the chiral order 
parameter scales as 

(W)(t,h) = h l l 5 F(z), (28) 

with z = th 1 ^ 6 , and the chiral susceptibility, defined as Xm(t, h) = d(ipip) / 'dm, 
via 

Xm(t,h) = - 6 0- l [F(z)-^F'(z)]. (29) 
The more familiar critical exponents S and (3 are related to y t and y% via 

P = $-ri and 5 = ^. (30) 

yt i - Vh 

The heights of the peaks of the susceptibilities scale themselves with the 
behavior 

Xr k ~™" 2m and Xt eak ~m- Zm (31) 

with z m = 2 - l/y h and z t = (y t - l)/yh + 1. 

The expected values for the critical exponents for the case of 0(4) sym- 
metry, 0(2) symmetry, and mean field exponents (MF) are listed in Table 
I, in the form of yh, yt, and the corresponding values of z m and z t . The 
0(2) symmetry exponents are also listed, because at finite lattice spacing, 
the exact chiral symmetry of the staggered fermion action is U(l) ~ 0(2). 
Only sufficiently close to the continuum limit does one expect to find 0(4) 
exponents. 
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Vh yt 


Z-m Zt 


0(4) 0.83 0.45 
0(2) 0.83 0.50 
MF 3/4 1/2 


0.79 0.34 
0.79 0.39 
2/3 1/3 



Table 1: Critical exponents for 0(4), 0(2) and mean field theory (MF). 
Taken from [|J. 

The calculated results for the exponents themselves, evaluated on different 
spatially sized lattices, are summarized in Table 2. Comparing Tables 1 and 
2, one sees that at this stage, no definitive statement about the symmetry of 
the underlying Lagrangian can be made from lattice gauge theory. This is 
an indicator that vital study in this field is still necessary to determine the 
underlying symmetry group conclusively. It is probably necessary to increase 
the lattice sizes and move to smaller masses. 





8 3 12 3 16 3 


Zm 

z t 


0.84(5) 1.06(7) 0.93(8) 
0.63(7) 0.94(12) 0.85(12) 



Table 2: Critical exponents, as a function of the lattice size. Taken from 
Ref.i. 



2.3.3 Meson screening masses and £/a(1) restoration 

One of the questions that has raised some theoretical interest in the last few 
years is whether the Ua(1) symmetry, i.e. the symmetry if) — > e* Q75 ?/>, which 
leads to the non-conserved current Aq that is given in Eq.(|2~0D is also restored 
at finite temperature, at some point. For three flavors, this occurs trivially. 
A demonstration of this, following Ref. p3| is given. 

In SU(3), the statement that Ua(1) is restored, implies that m n = m^. 
Since the masses of the particles are determined from the vacuum expectation 
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values of the appropriate meson-meson correlators, we need to show only that 



(<h(x)M<>)) = (0o(^)0o(O)), (32) 

where (j)^{x) = ip(x)i'j 5 \ 3 ip(x) is the correlator for the 7r and 4>q(x) = 
iji^Xoi/jfa) is that for the rj '. If one considers the specific axial transfor- 
mation 

ip(x) -> i/j'(x) = e il5{V ^ Xa - X:i) ^i/j(x), (33) 
then, after a little algebra, one finds that the composite fields transform as 




3 (x) - 0' 3 (x) = \I-Mx) + y 3^8- (34) 
The correlator composed of these composite fields itself then transforms as 
(M*)Mty) - (0' 3 (a;)0 , 3(O)) 



(35) 



The last two terms of this expression vanish, since the system is assumed 
to be S77y(3) symmetric. In addition, this implies that (0 3 (x)0 3 (O)) = 
(0 8 (x)0 8 (O)), so that Eq. fl3"5]) implies that 

(Mx)Mo)) = (<h(x)Mo)h ( 36 ) 

or that m v > = m^. 

In retrospect, it is simple to understand why the symmetry must be re- 
stored. Noting that mathematical constructions containing traces of fields 
preserve the symmetry, while determinants or antisymmetric functions vio- 
late it, one sees that the lowest order combination of fields that would violate 
Ua{1) would involve the completely antisymmetric tensor, and consequently 
contain three field combinations. Since one requires here only two field oper- 
ator combinations in order to construct a meson-meson correlator, this must 
be invariant in the chirally restored phase. This leads to the definitive 

statement: for n < Nf, all n-point functions in the chirally restored phase 
are Ua{1) invariant. 
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From the previous argument, it is evident that in SU(2) the situation is 
more complicated. There are two independent chiral multiplets in this case: 
(cr, 7?) and (7/, do). In lattice studies, the behavior of the masses of the ir and 
the do have been calculated. Here the integral over the correlators has been 
studied, 

XMi = [ d^iix^m) (37) 



for Mi = 7r or a , and the leading behavior of these correlators is assumed 
to be x^h ~ m ~M i - A pl°t °f the "screening masses" obtained in this fashion 
is shown in Fig. 2 as a function of 6/g 2 , with g the coupling in the QCD 
Lagrangian, which again represents increasing temperature over the region 
of the phase transition. It is interesting to note that the it and a have 
become degenerate: in this picture, this occurs at some temperature slightly 
larger than T c , with the a undershooting the it curve and approaching it 
from below. That the a meson undershoots the 7r curve is not expected from 
model calculations and may be a lattice artifact. This will be discussed in the 
following sections. One sees in Fig. 2 that the mass of the other scalar, the do, 
drops with temperature, but not as drastically as does the a. One observes 
that it does not become degenerate with the 7r and a over the temperature 
range indicated. Thus it does not appear from this particular calculation 
that Ua{^) symmetry is restored in SU(2). An alternate approach, however, 
using the scaling arguments of Brown and Rho EJJ], does however indicate 



a degeneracy at the transition temperature [4j. Thus, in this section once 
again, the question of the restoration of £7a(1) symmetry is not resolved. 



2.3.4 Bulk thermodynamic quantities. 

One of the most important contributions that lattice physics is able to provide 
are calculations of bulk thermodynamic quantities. In particular, the energy 
density and pressure densities are given by 

T 2 d 

e = ——\nZ (38) 
V dT K 1 

and 

p = T^\nZ (39) 
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Figure 2: Meson screening masses plotted as a function of (3 = Q/g 2 . The 
critical value /3 = /3 C is indicated by the vertical line. 



in terms of the partition function Z. In practice the pressure density 
is obtained from integration of the difference of the action densities at zero 
and finite temperature, 
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f/ 0o = N*J^d(3(S o -S T ). (40) 

Note that this quantity is defined in such a way that p/T^ =0 = 0, in contrast 
to setting the usual thermodynamic limit of Nernst, i.e. the entropy S(T = 



0) = |26| . While this does not affect anything that follows, one should bear 
this in mind when making model comparisons, as will be done in Section 3. 

In the following figures, we have chosen to illustrate the pressure and 
energy densities for lattice simulations that include quark degrees of freedom, 
rather than simply quenched QCD. In Fig. 3, we show the pressure density, 
plotted as a function of the scaled temperature, for four flavor QCD on a 
16 3 x 4 lattice. A comparison is made on varying the quark masses, and using 
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Figure 3: Pressure density, plotted as a function of the scaled temperature 
T/T c . (Taken from §.) 

quenched QCD, in the latter case with appropriate scaling of the number of 
degrees of freedom. One sees that there is a sharp rise in the pressure density 
at T = T c , and the curve tends to the Stefan-Boltzmann limit, but does 
not reach it over the temperature range (3.5T C ) shown. The deviation from 
the ideal gas limit appears to be too large to be described by perturbation 
theory, suggesting here that the perturbative regime occurs for temperatures 
T » T c . 



The energy density of four flavor QCD on a 16 3 x 4 lattice is shown in 
Fig. 4. In this case, the energy density remains close to the ideal gas limit at 
temperatures of the order of 3T C , but overshoots it and approaches it from 
above for finite values of the quark mass. Whether this is a lattice artifact 
or not is presently unclear. 
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Figure 4: Energy density plotted as a function of the scaled temperature, 
T/T c . (Taken from §.) 



In concluding this section, we see that lattice gauge simulations are reach- 
ing a point where one may obtain "exact" results that stem directly from the 
discretized QCD Lagrangian. These can be used as a guide for constructing 
simple models, and conversely, simple models and simple predictions based 
solely on symmetry considerations such as discussed heref] may be used as a 
guide for the interpretation of the numerical results. As has emerged here, 
there are still many questions that are open for study 

With this, we turn to a different approach which is regarded by its pro- 
tagonists as being an exact low energy representation of QCD, viz. chiral 
perturbation theory and investige what is known at finite temperature. 



3 Chiral random matrix theory |27| also falls into this category. 
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2.4 Chiral perturbation theory 



2.4.1 A brief introduction 

Chiral perturbation theory starts with the premise that an effective La- 
grangian for QCD at low temperatures can be written solely in terms of 
the observed baryonic (here mesonic) degrees of freedom, in such a way that 
global chiral symmetry is enforced. This is done in its most general form by 
collecting the mesonic degrees of freedom into the field 

U(x)=eM™ a T a /F), (41) 

where 7i a are the SU(2) pion fields, r a the Pauli matrices, and F the pion 
decay constant, and contructing a Lagrangian density that is ordered in mo- 
menta. Such an expansion for the Lagrangian only starts at 0(p 2 ), and 
must contain an even number of derivatives in order to be Lorentz invariant. 
Writing 

Cqcd - C eff = L% + C% + C% + . . . , (42) 
the lowest leading order term is 

4?/ = \Fhr(d,Ud^), (43) 

which, taken on its own, is the (non-renormalizable) sigma model. QCD, as 
we have already discussed, is however not completely invariant under chiral 
symmetry. There is an explicit breaking of the symmetry due to the presence 
of the current quark mass matrix. The symmetry breaking term is in general 
given as 

C sb = f(U,dU,...)xm°, (44) 

where m° is the (real and diagonal) current quark mass matrix. One incor- 
porates this into the effective Lagrangian by making not only an expansion 
in powers of the derivatives, but also in powers of m°, i.e. C s b ~ f{U) x m° 
to leading order. More precisely, this term takes the form (that is Lorentz 
and parity invariant) 

C sb = ^F 2 Btr(m°(U + U^), (45) 
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introducing the new constant B. This is generally included in the definition 



of Cf} f , i.e. 



df) = \F\i{dJJd»U ] ) + ]-F 2 Bti{m\U + U r )). (46) 
4 2 



In this reckoning, one can thus regard m° as being of 0(p 2 ). 

To make physical sense of the constant B, one may expand the field 
U = exp(i7r-r/F) in powers of the pion field ir. The symmetry breaking part 
of the Lagrangian then becomes 

C sb = K + <)B[F 2 - \* 2 + ^F~ 2 + ...} (47) 

The first term in this expansion gives the vacuum energy generated by the 
symmetry breaking. The second term generates the pion mass, while the 
higher order terms describe further interactions of the 7r fields. By direct 
analogy with the QCD Hamiltonian, we know that the derivative of Hqqd 
with respect to m' generates the operator uu. Thus the derivative of the 
vacuum energy with respect to the current quark mass gives the vacuum 
expectation value of this operator. Applying this to £ e ff, one has 

(0|uw|0) = (Q\dd\Q) = -F 2 B{1 + 0(m)}, (48) 

indicating that B is related to the condensate. Since the pion mass is given 
as 

m 2 n = (m° u + m° u )B{l + 0(m)}, (49) 
one obtains the Gell-Mann-Oakes-Renner (GOR) relation p8| , 

F 2 M 2 = (m° u + m° d )\(0\uu\0)\ (50) 

from Eqs.(^HJ) and (i9|), on eliminating B. 

To order p 4 , the effective Lagrangian would contain two additional inde- 
pendent terms in the event that no current quark mass were present, i.e. one 
would include two new terms 

Cf) f = h^tiid^d^U}) 2 + ^l 2 ti(d^d u U)tr(d fl U^d u U). (51) 

with new low energy constants l\ and Including the current quark mass 
matrix again to construct an explicit symmetry breaking terms requires the 



22 



inclusion of further additional terms, as was the case for C e jf. For most 
purposes, this is sufficient. However, to obtain the most general form from 
which all propagators can be derived, it is useful to introduce external fields 
into the Lagrange density. Here the essential additions are v^(x) and a^(x) 
that are vector and axial vector in nature and which can be regarded as being 
of order p. Then, using the original notation of Ref . [0J , the complete set of 
terms that contribute to dftt were worked out by these authors and found 
to be, for SU(3) 

L% = L 1 (V M f/ t V^) 2 + L 2 (V^ t V M f/)(V^ t V^) 

+ L 4 ( X U ] + X j U)(VU^U) + L & ( X U ] + X ] U) 2 

+ L 8 (xU j X U j + Ux ] U X + ) + L 10 (U^F^UF L , U ) 

+ H^FgFR^ + F£"F Lflv + H 2 ( x j x) + L 7 ( X U^ - U X ] ) 2 

+ L,(V,UV»U\ X U ] + U X ] )) + iLgiF^VVtf + F«VV U U) 

+ L 3 (V tl UVU^V v UV v .U^) (52) 

where using a different notation to Eq. fl5"T|) now, the low energy constants 
L\ to Liq, and Hi and H2 have been introduced. The angular brackets are 
a shorthand notation for the trace. In this expression, one notes that the 
covariant derivative that is constructed using the external field must now 
appear, 

V^ = <^-*K,f/}, (53) 
and F^ u is the field strength tensor constructed from the external field, i.e. 

Fjg L = ±d»a v =f d u a» -{[a", a"}. (54) 

Terms involving the current quark mass have been summarized into the field 
X = 2Bm, with m = (mJJ + Note that the low energy constants Lj 

become renormalized when physical quantities are calculated, as this theory 
is perturbatively renormalizable order by order. A certain number of such 
physical quantities that are measured in experiment must then be used to fit 
the renormalized parameters at a given mass scale. Given definite values for 
these constants, predictions of other quantities can then be made. 

Three ingredients are essential to any application that attempts to calcu- 
late quantities for chiral perturbation theory to a specific order. For example, 
should one wish to calculate to 0(p 4 ), the following steps must be taken: (1) 
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The general C e L of order p 2 is to be used at both the tree and one loop 

level. (2) The general ££jf of order is to be used only at tree level. (3) A 
renormalization program must be implemented to make physical predictions. 
The extension of this procedure to higher powers in p 2 is obvious. 

Let us look at a standard example for the derivation of the pion mass |29| . 
In what follows, we denote the low energy constants appropriate to SU(2) [] 
two flavors as being L[ 2 \ If one expands the Lagrangians CfL and dfjt in 
terms of the pion fields, one finds 

4// = jj[<9^7r-m 2 7^ 

(55) 

while 

- ^[32L« 2 » + IdL^-A ■ * + 0(0"). (56) 

The terms in C^jj that are of order 7r 4 contribute to physical quantities via 
one loop diagrams and one therefore does not need to consider these in a 
calculation to order 0(p A ). What is required however, are the one loop dia- 

(2) 

grams that are generated by £ e Jf- For a calculation of the the renormalized 
pion mass, however, one can avoid evaluating any diagrams at all by simply 
considering all possible contractions of two fields in these terms in C^)f, to 
arrive at an "effective" effective Lagrangian, that takes the form 

= I^vr • d„n - -m 2 n ■ n + )vr • vr 

eff 2 ^2 12m 2 V 7tJ 

+ Qpl&kdji - 5 u 5 k i)I(ml)(5 ij d> J 'n k d^TTi + bum^-K^j) 

+ \^,^Lf + 84 2) ] - \mln • ^[32^ + 16^]. 

(57) 

In obtaining this result, the Feynman propagator 

»A Fjfc (0) = ( 0|7r j (x)7r fe (x)|0) = S jk I(ml) (58) 

4 These can be simply related to the U of Eq. (|5l|), and the reader is referred to Q for 
explicit details. 
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has been introduced and is written in terms of the integral 

>K) = /"7 tSt^T = S-'W t2 ni - $)(ni)t- 1 (59) 



that is treated with dimensional regularization, d being an arbitrary dimen- 
sion. In addition, use has been made of the fact that derivatives of the 
Feynman propagator, defined as 

- d^d u iA Fjk (0) = (0\d f ,7r j (x)d L ,TT k (x)\0) = Sijl^(ml) (60) 

can be expressed in terms of the integral I{rn^) via 

WO - f j <gA*^ - fc^'M). («> 

Regrouping the kinetic and mass terms, Eq.([5TD becomes 

£*ft = ^-Ml + (164 2) +84 2) )||-3^/K)] 

- imV.7r[l + (32 4 2) + 164 2) )||-^/K)]. (62) 

By expanding this expression in powers of d — 4 and renormalizing the pion 
field as 7r r = Z~ 1//2 7r, with 



z. = i- ^(2Lf + 4 2) ) + Sf^t^i + -> ~ 1 " to 4 * + ln f 1- (63) 



one obtains the canonical form for the effective Lagrangian for pion fields, 

£eff = ^d^7l r d U 1T r - ^M 2 7T r • 7T r , (64) 

with the identification of the physical pion mass as 

M - =m2[i+ A 1,1 f - if ^ < 65 > 

Here L com (, = 2L^ r + L^ r — 4Lg 2 ' )r — 2L^ r . In the original paper of Gasser 
and Leutwyler M , M% was not obtained in this fashion, but rather from the 
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expansion of the Fourier transform of the axial vector correlator, which has 
the form 



J%(P) = ij dV*-» ) (0|7Vt;M^(y)|0) 

= ^Jif$ + --^ < 66 > 

where A 1 (x) = ip^x^^^ip^x). From this expression, the corresponding 
expansion for F v has also been obtained. 

2.4.2 Cool chiral perturbation theory 

The evaluation of the condensate density at finite temperature was first car- 
ried out by Gerber and Leutwyler ||. In their calculation, which involves 
C^ff and £e//> they find that the first term in the behavior of the condensate 
with temperature is quadratically decreasing, i.e. 

(55) = (0|5,|0> T .o[l - ^- 3^1-5^'^ + 0(7*)]. (67) 

This is a result that has been obtained under the assumption that quarks 
are massless, i.e. in the chiral limit. A q is a scale factor constructed from 
the renomalized low energy constants, and is expected to be of the order of 
A q = 360. .580 MeV. 

In a recent publication, Toublan || has investigated pion static properties 
with the aim of obtaining 0(p 6 ) accuracy in all quantities and to then verify 

the Gell-Mann-Oakes-Renner (GOR) |28[ relation at finite temperature. To 

(2) 

do so, the tree, one loop and two loop diagrams of C e L are required, the tree 
and one loop graphs of C^h are required plus the tree level graphs of jQ//- 



In doing so, the result of Eq.(67) has been reconfirmed. In addition, the mass 
M n (T) and pion decay constant as a function of temperature are also evalu- 
ated, using the finite temperature axial vector correlator. In total, thirty-six 
Feynman graphs contribute to the correlator at this order! However, in the 
chiral limit, one is still lucky enough to have simple analytic forms for the 
temperature dependence. One finds 

M ' (T) -1+ T ' ^ ln AM + OfT 6 ) (68) 
~M~ - 1 + 24i^ _ 36^ ln T + 0iT 1 (68) 
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while 



and 



Re[Fl(T)f. T 2 T 4 A 



i-T^ + ^i^ + ^T 6 ). (69) 



^^^^M.T.a^), (70) 

where A^t.a are various scales, whose sizes are determined by the renormal- 
ized couplings L\. . . L[ Q that are a function of scale. They are determined 
numerically to be A M ^ 1.9 GeV, A T ~ 2.3 GeV, and A A = 1.8 GeV. Note 
that, at finite temperature, there is a separation of "temporal" and "spa- 
tial" pion decay constants. This comes about since Lorentz invariance is not 
maintained in a heat bath and the the singular part of the axial two point 
function takes the form 

where 

f (q,T) = q F t (q,T) f t (q,T) = q t F s (q,T), (72) 

with i = 1..3, and the decay constants F^ are defined as 

Ff(T) = F s>t (q,T)\ ?=0 . (73) 

The GOR relation is modified so as to read |J 

Um MllT)MFl(T)Y = _ t 
m-^o m(qq)T 

For this reason, we show graphs for M 2 (T)/M 2 and Re[F^(T)] 2 /F^, as a 
function of temperature in Figs. 5 and 6. The tree level result is given (dotted 
curve), together with the one loop computation (upper dashed line in Fig. 5, 
lower dashed line in Fig. 6) and the two loop approximation (solid curve). In 
both of these figures, a non-zero value of the quark mass has been assumed 
for these curves. In the chiral limit, one finds the lower (upper) dashed 
curve in Fig. 5 (6). What is evident from these two figures, is that chiral 
perturbation theory is not converging and appears to provide an oscillating 
series for these quantities. Thus for larger values of the temperature, T > 100 
MeV say, one sees that the pion mass decreases with temperature in the two 
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Figure 5: The pion mass, scaled by its value at T = is shown as a function 
of the temperature. 

loop approximation, in contradistinction with the one loop result, the lattice 
results of the last section, and also in contradistinction with the model results 
obtained in the Nambu-Jona-Lasinio model, which will be presented later on 
in the following section. Convergence at temperatures in this range appears 
to be problematic, which is perhaps an indication that the series is at best 
asymptotic, or changes its nature due to the onset of the phase transition. 
In this range, one expects non-analytic behavior and it is unreasonable to 
expect a perturbation analysis to succeed. These curves clearly indicate that 
ChPT at finite temperature can at best be regarded as cool, so that the 
fundamental behavior at low temperatures sets a constraint on the finite 
temperature behavior of would-be effective models. 
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Figure 6: The pion decay constant, scaled by its value at T = is shown as 
a function of the temperature. 

3 The Nambu— Jona-Lasinio model 

The Nambu- Jona-Lasinio (NJL) model has been reviewed in detail by several 
authors from different viewpoints || P, |I0||, and consequently I do not wish to 
present any detail of this model other that a basic introduction here. Rather 
the purpose of this chapter is to illustrate that with the simple equations 
requiring little computational time, one can reproduce all the main features 
of the static properties that have been so arduously extracted from years 
of labor on the lattice. It is extremely encouraging to have a simple model 
that can be handled semi-analytically - one gains a tremendous amount of 
insight into the actual functioning of the mechanism of dynamical symmetry 
breaking and the consequences thereof. 

Nevertheless, the NJL model is simply a model - in contrast to the results 
of the previous section, which are regarded as "factual" , this section can only 
give model-dependent results. Accordingly it is only equitable to indicate, in 
addition to the successes provided by this approach, the failings also. These 
become obvious when examining bulk thermodynamic properties, such as 
pressure, energy and entropy densities, and will be discussed in what follows. 
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We shall then turn to dynamic properties, and examine the temperature 
dependence of scattering amplitudes in the quark-antiquark channel, which 
displays a divergence which we term critical scattering, in analogy to the 
phenomenon of critical opalescence that is observed in light scattering. 



3.1 Order parameter 

We first consider the order parameter for the chiral transition that is obtained 
from the NJL Lagrangian, which, for two flavors of quarks, is taken to be 

Cnjl = 1>(x)(i - mMx) + G[{^f + (^ 5 rtfj)% (75) 

where G is a dimensionful coupling strength, and m denotes the common 
current quark mass for u and d quarks. For three flavors of quarks, we use 

a=0 

- K{det^(l + 75)^ + det^(l-7 5 )7/;}. (76) 

Here G and K both are dimensionful coupling strengths and = diag(mQ, m[J, 
The self-energy, in the mean field approximation, that corresponds to the low- 
est order term in an expansion in the inverse number of colors iV c |I2 , |l3|j , is 
given as[] 

m = m -2G((^», (77) 
and the condensate is given explicitly as 

N Nr r A r? 

m)) = ^- I l ^r[l-r(p,»)-f + (p,»)}, (78) 

7T Z JO tip 

with 

/ ± fe' u ) = T7T L F z Vv ^ 

[1 + exp /3{E P ± fi)\ 

One sees that the condensate is directly proportional to the value of the 
dynamically generated mass, in the event that the current quark mass is zero. 



5 Since the coupling strengths turn out to be large, GA 2 ~ 2, an expansion in the 
number of couplings is inadmissable and an alternative expansion scheme must be used. 
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Although the situation is more complicated in SU(3), where the dynamically 
generated quark masses satisfy coupled equations, 

o GN c A , \ KN 2 C 
mi = m i —miA(m u p-) H —mjm k A{mj, \ij)A{m k , fjL k ), 

i^j^k (80) 

and the function 

A(m f ,fi f ) = ^Ve 4 ^ / \ 9 ^ 2 (81) 

V ' /7 /? V ibKA (27r) 3 (^ n + /i / ) 2 -E2 v ' 

is proportional to the condensate density for a specific flavor, 

A(m f ,iif)~({W,))f, (82) 

the dynamically generated quark masses are equivalently order parameters of 
the phase transition, and we therefore plot these. They are shown here only 
for the SU(3) case, in Fig. 7, for a finite value of the current quark mass |30| . 
As expected, the phase transition that occurs in the chiral limit is washed 
out and becomes a cross over. Another feature that emerges in this model 
is that the strange quark mass remains large, even at temperatures T ~ 300 
MeV, and does not reach its current mass value of 150 MeV until T » T c . 



3.2 Meson masses 

The meson masses for the scalar and pseudosalar sectors are determined via 
the well-known method of evaluating the quark- ant iquark scattering ampli- 
tude in the random phase approximation, and searching for poles of this 
function. This involves knowing only the irreducible polarization function 
that one can construct from a single quark loop, the details of which can be 
found, for example, in || |TI| |51J. One finds the masses that are shown in 
Figs. 8 and 9 for the pseudoscalar and scalar sectors, respectively. In Fig. 8, 
2m q is plotted in addition to m n . The point at which these two curves cross 
is called the Mott temperature, For T > Tm„, the pion is no longer a 

bound state, but is a resonance, with a finite width that is not shown here. 
Similarly we have plotted m q + m s , from which the kaonic Mott temperature 
T Mk is defined. For T » T Mk , the kaon is also a resonance with a finite 
width 
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Figure 7: Temperature dependence of the constituent quark masses. The 
solid line refers to the light quarks, the dashed to the strange quark [RO . 



These graphs deserve some comment. Firstly let us compare them with 
the figure showing the meson masses obtained via lattice gauge theory, Fig. 2. 
We note first that there are some fundamental differences in obtaining these 
graphs: (a) Figs. 8 and 9 show so-called pole masses, while Fig. 2 gives 
screening masses. Nevertheless, it has been shown that, in the NJL model, 
the temperature behavior of screening masses and pole masses is qualitatively 
similar [31], although quantitatively somewhat different. Since we cannot 
hope for any quantitative agreement at this stage, it is justifiable to make a 
comparison, (b) The NJL model calculation shown is for SU(3), while the 
lattice calculation is SU(2). With these points in mind, note that the a and 
7r mesons from figures 8 and 9 become degenerate at high temperatures, as 
observed also in Fig. 2. However, there is no undershooting of the a meson. 
The meson labelled of Fig. 8 corresponds to the a$ of Fig. 2. Here 
one observes qualitatively the same behavior, i.e. that this scalar meson 
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Figure 8: Temperature dependence of the pseudoscalar meson masses, as well 
as that of 2m q and m q + m s . |30] . 



1000 



> 



500 




100 



200 
T (MeV) 



300 



Figure 9: Temperature dependence of the scalar meson masses and 2m q |30 
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also decreases strongly in the phase transition region. Thus one has an 
aesthetically pleasing agreement between the NJL model and the results 
obtained by lattice gauge theory for meson masses at this level. 

A direct comparison of Fig. 7 with the results of chiral perturbation the- 
ory, i.e. with Fig. 5 is problematic. We simply make some comments: (a) the 
physics underlying Figs. 5 and 7 is completely different. Fig. 5 is obtained 
by constructing meson loops (the mesons are regarded as structureless point- 
like objects), while in Fig. 7, the pion is constructed from a quark- ant iquark 
loop. Meson loops form corrections to this calculation and would be of the 
next order in 1/N C . Such corrections have in fact been evaluated, and it 
has been found that the leading order T 2 dependence of Eq. (|68|) is recovered 



32| . The fact that in the final analysis the curve of M w as a function of 
temperature is finally decreasing for chiral perturbation theory in Fig. 5, is 
in strong contradiction to both Figs. 2 and 8. 



3.3 Bulk thermodynamic quantities 

In the last subsection, we have indicated the successes of the NJL model in 
calculating the order parameter and masses as a function of temperature. In 
this subsection, we turn to bulk thermodynamic quantities. Here we will see 
that the model does not do as well, and that the lack of confinement makes 
itself strongly evident at low temperatures, while the cutoff of the model 
is a hindrance at high temperatures. We start with the thermodynamical 
potential Q, calculated in the grand canonical ensemble. Given an interaction 
between fermions that is 4-point in nature such as in Eq.([75|), Q can be 



calculated quite generally as [10, 33 



r 1 d\ 1 r rl 3 n 1 

n = n ° + Jo T 2 J ^pT,^P^nVm^>n,P)S X (u n ,p)}, (83) 

where Qq is the thermodynamic potential in the absence of interactions, 
and Y, X an d S x designate the Matsubara self-energy and Green function 
associated with the system. The superscript A refers to the fact that both S 
and £ are to be evaluated with the introduction of an artifical coupling that 
multiplies the interaction Lagrangian £j nt . The Matsubara frequencies for 
fermions circ, £is required, odd, i.e. v n = (2n + 1)tt//3, with n — 0, ±1, ±2 ± 
3.... 
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For the NJL Lagrangian of Eq. (|75|) , in the mean field approximation, it 
is not necessary to apply Eq.(^). A straightforward calculation gives 

n = n, = -2n c n, J ^E r - 2 -Mi j ^ M ne-^\i +e -^-»X 

(84) 

with E£ = p 1 + m 2 . As is evident from the label q, this appears to be a 
thermodynamic potential generated solely by the quark degrees of freedom. 

We note also that the thermodyamical properties can only be measured 
relative to the physical vacuum, 

^» = fi(T = 0,/i = 0,m(0,0)), (85) 

which, for the mean field approximation, corresponds to 

(«*W=^^-2^/^. (86) 

To introduce mesonic degrees of freedom, it is necessary to go beyond 
the mean field approximation to include the next set of terms in the 1/N C 
expansion. The self-energy in this case includes effective interactions in both 
the scalar and pseudoscalar channels |34| , 

1 r d 3 q 

+ ij 5 TS(v n f, q)i^rV^{y n -v n ,,p-q)], (87) 



and is constructed on summing the Fock and infinite RPA series that con- 
tribute to the self-energy in this order. Here 

V&(w, q) = -2GX[1 - 2GXU M (u, q)Y\ (88) 

and the irreducible polarization in the mesonic channel 

Tl M (uJ m ,p) = [ -^-^^^TrT M S(u n + LO m ,q + p)T M S(u n , q), (89) 

is determined by the vertex Tm for that channel. Inserting Eq.(^) into 
Eq.(^) yields the fluctuating part of the thermodynamic potential, 

fi /< = E^f / (^s^E^'Mi -2Gn M K,p)]. (90) 
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The nature of this term is revealed on performing the frequency sum. One 
has 

n fl = n v + n ff , (91) 

where, for each species M = ir or a, 



d i-2Gn M (o; + «6,p) 

du l-2GU M (oo-ie,p}' v ; 

Some analysis shows that a simple approximation for the polarization near 
the pole, i.e. 

1 - 2GU M (oj,p) = (lu 2 - E 2 M ) x const (93) 

leads to 

= N M J ^[\e m + i ln(l - (94) 

exactly as one would expect for the thermodynamic potential given bosonic 
degrees of freedom. The pole approximation is however insufficient, as one 
integrates over all energies, and in practice, the fact that the bound states 
also become delocalized resonances at the Mott point must also be accounted 
for. This has been done in introducing phase shifts in each channel ||34|| . 

In order to calculate the pressure, the physical vacuum given by Eq.(p5|) 
must be reevaluated to include a term from Qfi. One now has 

^tr = + (nftV ( 95 ) 

and the pressure density is now 

p = -n q - n fl + nCr- (96) 

In Figs. 10 and 11, the pressure and associated energy densities evaluated 
from this thermodynamic potential are shown. 



In Fig. 10, one sees that the lower curve, corresponding to the mean field 
approximation calculation represents the quark degrees of freedom. There is 
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Figure 10: The pressure density, scaled by T 4 , is shown as a function of tem- 
perature. The lower curve is for the mean field case only, the upper includes 
fluctuations (mesons). The vertical line indicates the critical temperature 
and the horizontal one the Stefan-Boltzmann limit for an ideal quark gas 
PI. 



an appreciable pressure that arises from this term, i. e. from the quark degrees 
of freedom, for temperatures T < T c , which is indicated by the vertical 
line. Including mesonic degrees of freedom rectifies the behavior at small 
temperatures, but still leaves a large intermediate range of temperatures 
T < T c that is dominated by these unphysical quark degrees of freedom. This 
is thus a direct consequence of the missing feature of confinement. The sharp 
rise in the pressure density shown in Fig. 3 cannot be modelled by a non- 
confining theory. At high temperatures, T > T c , there is a small contribution 
from the mesonic degrees of freedom, that exist as correlated states with a 
finite width in the plasma. The main contribution arises however here from 
the quark degrees of freedom. The actual value obtained for the pressure 
density underestimates the Stefan-Boltzmann limit (shown as a horizontal 
line), since there is a finite cutoff on the quark momenta. Relaxing this 
constraint would lead to the pressure density approaching a constant. 



Similar comments can be made for the energy density: the intermedi- 
ate temperature range 50MeV< T < T c is dominated by quark degrees of 
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Figure 11: The energy density, scaled by T 4 , is shown as a function of tem- 
perature. The lower curve is for the mean field case only, the upper includes 
fluctuations (mesons). The vertical line indicates the critical temperature 
and the horizontal one the Stefan-Boltzmann limit for an ideal quark gas 
0- 

freedom, indicating the lack of confinement. The high temperature values 
T > T c do not approach the Stefan-Boltzmann limit, due to the cutoff. 

In concluding this subsection, one sees that one needs to include confine- 
ment in some fashion in order to be able to regain the lattice picture. From a 
thermodynamic point of view, the high temperature regime about T ~ T c is 
probably best described by the model, in the sense that only quark degrees 
of freedom plus correlated mesonic states are present. In the next subsec- 
tion, we thus study elastic quark-antiquark scattering about this point and 
indicate that a divergence occurs in the cross-section at T = Tm and that 
the phenomenon of critical scattering as a consequence of the chiral phase 
transition is observed. 

3.4 Critical opalescence in the quark-antiquark chan- 
nel 

In this section, we examine the behavior of the quark-antiquark scattering 
amplitude in the NJL model in the vicinity of the Mott temperature, which 
replaces the critical temperature when finite current quark masses are used. 
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In SU(3), there are seven independent processes out of a total of fifteen for 
quark- ant iquark scattering, taking isospin and charge conjugation symmetry 
into account. These are listed in Table 3. Mesons that can be exchanged in 
the s and t channels, as are given by the Feynman diagrams of Fig. 12 are 
also listed. 



Process Exchanged mesons (s channel) Exchanged mesons (t channel) 

ud — > ud it, <j n 7i, r\, rj', ov, a, a' 

us — > us K, <7k Vi rfi 'i °' 

uu — > uu it, rj, rj' , a n , a, a' it, rj, rj' , a n , a, a' 

uu — > dd 7r, rj, rj' , a n , a, a' it, a n 

uu — > ss rj, rj' , a, a' K, a K 

ss — > uu rj, rj' , a, a' K, a K 

ss — > ss rj, rj' , a, a' K, a K 



Table 3: Independent processes for qq scattering. 



The transition amplitudes can be written as 

-iM s = 5 CUC2 5 C3:C4 v(p 2 )Tu(p 1 )[iV^(p 1 +p 2 )}u(p 3 )Tv(p 4 ) 

+S cl ,c 2 S C3 ,c 4 v(p2)(il5T)u(p 1 )[iV^(p 1 + p 2 )}u(p3)(i^T)v(p 4 ), 

(97) 

and 

-%M t = 5 CuC3 5 C2:C4 u(p 3 )Tu(p 1 )[W^(p 1 -p 3 )]v(p 2 )Tv(p 4 ) 

^c 1 ,c 3 Sc 2 ,c 4 u(P3)(il5T)u(p 1 )\iV[(p l - P3)}v(p 2 )(il5T)v(p A ), 

(98) 

S P 

where T selects the isospin eigenvalue for a particular channel, and D s ^ is the 
s or t channel, scalar or pseudoscalar quark- ant iquark scattering amplitude, 
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Figure 12: Feynman diagrams for elastic qq scattering within the NJ1 model 
in an expansion to lowest order in 1/N C . 



and which can be constructed from the corresponding polarization function. 
It has a simple form, for example || 

where G e ^ is an effective SU(3) coupling strength in the pionic channel |§. 
The differential cross section is constructed in the usual fashion as 

dt = 16n[s - (m u + m s ) 2 }[s - (m u - m s ) 2 ] 4iV c 2 ^ ^ Ma " Mt ^ (1 ° 0) 

while the total cross section is evaluated as 

da 



a 



J dt ^[1 - f F (PE 3 )][l - MPEJ], (101) 



introducing a Fermi blocking factor for the final states. Here Ef = pf + mf, 
where i = 3,4. 

In Fig. 13, we show the total cross section for light quarks in the ini- 
tial state, as a function of y/s, at a temperature T = 215 MeV, which lies 



40 



45 - 

40 - 

35 - 

^30 - 

JEU5 - 

b 20 - 

15 - 

10 - 

5 - 

- 



~i 1 1 1 1 1 1 1 1 — 

rid —> ud 

uu — > uu 

mi —¥ dd 

uu — >• ss 

T = 215 MeV 



_l I I u 



0.2 0.4 0.6 0.8 1 

yfs (GeV) 



1.2 1.4 



Figure 13: Total cross section for qq scattering with only light quarks in the 
initial state, shown as a function of y/s, at T = 215 MeV. 



slightly higher than the pion and kaonic Mott temperatures, = 212MeV 
and Tm k = 210MeV. Both pions and kaons are sharp resonances now. At 
higher values of the temperature, these become broader resonances in the 
cross-section, as shown in Fig. 14. At the Mott temperature itself, when 
quarks bind into hadrons, the intermediate states in the s channel give rise 
to infinite cross sections at threshold. This feature, which also appears in 
other processes like tttt — > Tin ||35|| , 7T7 — > 7T7 |]36| or qq — * 77 ]37 , [3~8| is akin 
to the phenomenon of critical opalescence. This has been discussed in some 
detail in Ref.[l39|l, and the interested reader is referred to this. 
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Figure 14: Total cross section for qq scattering with only light quarks in the 
initial state, shown as a function of y/s, at T = 250 MeV. 



4 Non-equilibrium formulation and transport 
equation 

The considerations of the first two sections discussed properties of chiral 
systems in equilibrium. If it were possible to measure any of the associated 
changes at the phase transition temperature, there would be no need for 
further discussion. However, because of the nature of confinement, we are 
unable to observe critical scattering directly, nor any of the other dramatic 
changes in pion properties. One tool for examining quark matter is via heavy 
ion collisions, and as such, over the short time scales over which collisions 
occur, it is unclear whether both thermal and chemical equilibrium can be 
reached during a collision. For this reason, we wish to investigate what the 
effects are of a condensate that changes with the medium, as well as medium 
dependent cross sections in a non equilibrium scenario. 

There are several formal aspects that have to be understood before one 
can attempt actual collision simulations. Firstly one can set up an exact 
formal description of a relativistic fermionic system that is out of equilibrium 
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via the method of Schwinger and Keldysh. From a heuristic point of view, 
however, we have a good understanding of the classical Boltzmann equation, 
so that it is important to establish a link between the two from which one can 
then go further. In doing so, one generally has a field theory with retarded 
and advanced Green functions. However, if we examine the collision term of 
the Boltzmann equation, we see that we require cross sections. However, we 
only know how to calculate these using causal Green functions. So we have 
to find a link telling us which level of approximation requires which Feynman 
graphs. 

The content of this lecture is summarized briefly in the next paragraphs, 
(a) We wish to start from a chosen Lagrangian that gives a microscopically 
correct description of the world, and to formulate a non-equilibrium theory 
via a matrix of Green functions S^Xi, £2) (* and j will be defined later!). 
This matrix of Green functions satisfies a matrix form of the Schwinger-Dyson 
equations, which as usual, cannot be solved exactly, (b) Some technical aid is 
required at this point. A centre of mass variable X = (xi+Xz) /2 and relative 
coordinate u = X\—X2 are introduced, and one Wigner transforms the matrix 
of Green functions. This is simply a Fourier transform with respect to the 
relative coordinate u. At this point, the equations are still exact, (c) Now 
one seeks methods of solution. For a fermionic system, the exact method 
would involve making a spinor decomposition of the Green functions, and 
we would have 32 coupled equations to solve! This is simply too difficult, in 
particular for an expanding system, for which spatial gradients are important, 
and so we turn rather to making the quasiparticle assumption, which, coupled 
with an expansion in powers of h, leads to the well-known kinetic theory of 
Boltzmann, here in relativistic form. 

All that has been discussed is quite general for any fermionic theory. 
Using the Lagrange density of the Nambu-Jona-Lasinio Lagrangian with an 
expansion in 1/N C illustrates how extensions to the standard binary collision 
forms in the Boltzmann equation come about, and clears the issue of the 
content of Feynman graphs for the cross sections that occur in the Boltzmann 
equation. 



4.1 Closed time path — Schwinger-Keldysh formalism 

There are several excellent texts that exist that cover the basics of the 



Schwinger-Keldysh formalism jll], [LLj for Green functions not in equilibrium. 
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Detailed reviews using path integrals can be found in H2J, while the more 



standard operator approach is to be found in |44], |45[ . Most confusing in 
this subject is simply notation: All the listed references use different ones. I 
shall conform to that of Landau^, which is particularly transparent in setting 
up rules for a perturbative diagrammatic expansion. 

Central to the problem of non-equilibrium systems is that the description 
via a single causal Green function alone, is inadequate. One requires the four 
Green functions, 

ihS c (x,y) = (Tip(x)$(y)} = iHS~~(x,y) 

ihS a (x,y) = (filj(x)ij(y)) = ihS ++ (x,y) 

ihS > (x,y) = (%l)(x)4)(y)) = ihS + ~(x,y) 

ihS<(x,y) = -(ij(y)ij(x)) = ihS- + (x,y), (102) 

i.e. the causal and acausal propagators S c and S a , S > and S < . In Eq. ( |102j ), 
T is the standard time ordering operator, 

T(0(x)0(y)) = 6(x - y )O(x)O(y) - 6(y - x )O(y)O(x), (103) 

and T the anytime ordering operator, 

f(0(x)0(y)) = % - x )O(x)O(y) - 6(x - y )O(y)O(x). (104) 

On the right hand side of Eq.( |102| ), the superscripts ij = +, — have been 
introduced (these were mentioned in the introduction to this section). This 
is an arbitrary but useful convention for constructing a matrix notation for 
summarizing the Green functions, 

It is automatically achieved by introducing the closed time path of Fig. 15, 
and setting the fields that occur in the Green function on the ith or jth 
branch respectively. 



6 This differs from the labelling of by a minus sign. Off-diagonal self-energies also 
differ by a minus sign. 
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Figure 15: Closed time path on which the Green functions are defined 



There are many interlinking relationships that follow simply from the 
definition of the Green functions. For example, S c and S a are related to S > 
and S < via 

S~~(x,y) = 9(x -y )S + ~(x,y) + 9(y -x )S~ + (x,y) 
S ++ (x,y) = 9(y -x )S + -(x,y)+9(x -y )S- + (x,y). (106) 

All four Green functions are not independent, since 

S~~(x : y) + S ++ (x,y) = S~ + (x, y) + S + ~(x, y) 

= S K (x,y) = ~(mx)My)]), (107) 

defining the Keldysh Green function. In addition, one can define the retarded 
and advanced Green functions 

ihS R (x,y) = 9{x G -y Q )({^{x)3{y)}) 

ihS A (x,y) = -e(yo-x )(Mx),$(y)}), (108) 

which are also related to the S lj via 

S R (x,y) = S~~(x,y) - S~ + (x,y) = S + ~(x,y) - S ++ (x,y) 
S A (x,y) = S-(x,y)-S-+(x,y) = S- + (x,y)-S+ + (x,y), (109) 

which can also be verified directly from the definitions of these functions. 
One could consider working with the matrix of independent functions 

£'=(c°R £k), (HO) 

but I will not do so in this chapter. Nevertheless, the retarded and advanced 
Green functions play a special role. Due to their simple analytic structure, 
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Figure 16: Dyson equation for one component of the matrix of Green func- 
tions. 



plus the fact that the equations of motion that they satisfy (see Eq. ( |117| ) 



later!) are closed, means that one usually can find a simple analytic form for 
these functions. 

The matrix of self-energies is defined now via the Dyson equation, 



S_(x,y) = Sl_(x,i/)+ I d 4 zd l wS (x, w)T,(w, z)S(z, y) 

= S^_(x,y) + J d 4 zd 4 wS{x,w)Uw,z)S^{z,y). (Ill) 

Pictorially, one can for example examine one element of this equation - say 
S ++ . The equation that this function satisfies is given in Fig. 16, using an 
obvious notation. Thus one sees that all components of the self-energy are 
in fact required in order to evaluate one single component of S_. 



From the Dyson equation, one can derive the equations of motion for the 
components of S_, which are summarized as 



(ih ft x -m )S(x,y) =qj> (x - y) + / d zaj:(x, z)S(z, y), (112) 

where 

a-(J-i)- < 113 ) 
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By defining the retarded and advanced self-energies as 

Z R = £— + £-+ 

S A = IT" + £+-, (114) 
one finds the corresponding Dyson equations for 

SL(x,y) = & { °\x,y) + J d 4 zd 4 wS!_ {0) (x,w)Yi(w,z)SL(z,y), (115) 

from which one sees that the Dyson equations for S R and S A are individually 
closed, 

Sff(x,y) = Sff°\x,y) + J d 4 zd*wSff {0) (x, w)tff(w, z)S^(z, y) 

(116) 

with corresponding equation of motion 



(ih <ft x - m ) a/3 Sf^ A (x,y) = 5 a ~ / 5 4 (x-y) + J d A zY%f{x, z)Sf^ A (z,y), (117) 

while the equation for the Keldysh function is integrodifferential, 

(ih @ x - mo) S K = J (H K S A + Y, R S K )d 4 z. (118) 

For a free particle, it is useful to note that the solution for the retarded and 
advanced functions follows immediately as 

S R > A (P) = 2 * + 2 ™°. ■ (H9) 

4.2 Transport and constraint equations 

Of the matrix of Green functions, consider only the equation of motion for 
S~ + (x, y) that follows from Eq.( |112| ). This is 



(ih ft x -m ) a pSp+(x,y) = j d 4 z (x, z)S p +(z, y) + Zj(x, z)S^+(z, y)) 

(120) 

In a similar fashion, one can derive the equation of motion 

S^fav) (-i0y - W)^ = J d 4 z{-S-p(x,z)J: p + (z,y) 

-S-J(x,z)^( Z) y)}. (121) 

47 



It turns out to be slightly more convenient to cast these equations in an 
alternative form, using the relations Eq. (|109|) between the Green functions, 
and similar ones for the self-energies. We write 

(ih @ x - m ) Q/3 y) = J d 4 z{S;+(x, z)S^~(z, y) - z)S^(z, y) 

+E* z)Sp + (z, y) - S-+(x, z)S*(z, y)} 

(122) 

and 

Srf{x, V) (~i0y - mo)^ = J d 4 z{ - S^ix, z)Y,~p + {z, y) 

+ S-+(x,z)Z^(z,y)}, (123) 

It is now a tedious technical task to Wigner transform Eqs.( |122| ) and ( |123| ). 
We illustrate this on a simple example and then simply give the final result. 
Introducing relative and centre of mass variables u = x — y and X = ^(x + y), 
the Wigner transform of S(x, y) is defined to be 

S(X,p) = J d\e^ h S(X + \,X- |). (124) 

To Wigner transform say the first term on the left hand side of Eq. 
requires an integral of the form 

d 4 ue ipu/h d^f{x,y) 

Id d s 1 „ 1 



/ d A ue ipu '\-^- - ^-)f(X + -u,X- -u) 
J 2dX^ du^ JJK 2 ' 2 ' 

I d 4 ue ipu/h f(X + -u,X- -u) 
2 dX„ J M 2 ' 2 J 



^f + j)/(A». (125) 



Similarly one can show that 



dZf(x,y) (-i^ + \d»)f(X,p) (126) 
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f(y)g(x,y) - / (X) exp [ I g(x, P ) (127) 

/(■'•)//(.r.//) - /(A")oxi> | | //(X./.) 1128) 

d 4 zf{x,z)g{z,y) -> /(A~,p) exp I -^( / - / 9 )\g(X,p) 



(129) 



need be made on Wigner transforming the product functions on the left hand 
side of the last equations. Applying these relations to Eqs.( |122j ) and (|123|) 
leads to the rather complex forms for the equations of motion, 

^{—fr - ^)--o}S- + (A>) = 

Z- + (X,p)AS + -(X,p) - Z + -(X,p)AS~ + (X,p) 
+ Z A (X,p)AS~ + (X,p) - X-+(X,p)AS R (X,p) 

(130) 

and 

S-+{X,p){-ifyf^^ + % ^) -m } = - S R (X,p)A^(X,p) 

+ s- + (x, P )Ai: A (x,p), 

(131) 

A = exp(-^(-^^-^-^)V (132) 

Now subtracting and adding these resulting equations, one arrives at two 
futher equations, which we identify as the transport and constraint equations 
respectively: 

V {7 "'^ }+[i, '^ +(x ' p)]=/ - (133) 

and 

T [7M '9A^ ] + { ^ m °'^ +} = /+ - (134) 
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In these equations, the terms that occur on the right hand side are decom- 
posed into three types of contribution, one containing at least one retarded 
function, one with at least one advanced function and a further term with 
neither, which in the semi-classical limit is the origin of the collision integral. 
Explicitly, one has 

I T = hoii (135) 



with 



and 



T 1 T ' 

I coll = E- + (X }P )AS + -(X,p)-E + -(X lP )AS- + (X,p) 

= 4 g ot-4X (136) 
/? = -Z-+(X,p)AS R (X,p) ± S R (A»A£- + (A» (137) 

I* = Z A (X,p)AS- + (X,p) T 5- + (X,p)AS A (X,p). (138) 

Equations (|133| ) and (|134j) are the central, exact equations that describe the 
non-equilibrium evolution of a system of interacting quarks. To actually see 
that these are in fact transport and constrint equations known from Vlasov of 
Boltmann theory requires some (hard) work. This follows only under certain 
approximations, and of course one needs some model in order to specify the 
interactions. For this purpose, we will use the Nambu-Jona-Lasinio model. 
Before doing this however, note that an exact solution of Eqs. (|133|) and (|134 ) 
follows formally on making a spinor decomposition, 

- ihS~ + — F + i l5 P + y% + 7 M 75^ = ^ V (139) 

The equations for the projected functions F ~ tiS~ + , P ~ trj 5 S~ + , . . . , 
form a set of 16 times 2 coupled equations that need to be solved simulta- 
neously. This is not only a formidable task from the computational point of 
view, it also offers at present little physical insight. 

For reasons of simplicity, therefore, we introduce the quasiparticle ansatz 
that contains the quark and antiquark distribution functions f q (X,p) and 
fq(X,p), and which puts these on their mass shell, 

S- + (X,p) = 2ni^±^[5( Po - E p )f q (X,p) - 5(p + E p )f q (X, -p)] (140) 

with f q} q = 1 — f q> g. Similar expressions can also be easily written down for 
the remaining components of the matrix . 
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4.3 The Vlasov equation for the NJL model 



At this point, one cannot go futher unless one specifies a theory or model from 
which the self-energy can be calculated. A four point interaction like that of 
the SU(2) NJL model is particularly simple to handle because the Feynman 
rules are particularly simple: (a) a directed line represents a fermion. The 
signs attributed to the beginning (i) and end (J) of the line reflect in the 
Green function iS^ 1 to be associated with the line, (b) an interaction line 
can have only a single sign on both of its ends. If the sign is ±, it is to be 
translated as ±iV, with V being the interaction strength. In the NJL model, 
this is V = -2G. 

According to these rules, in the Hartree approximation, it follows immedi- 
ately that £ + ~ = £ h = 0, so that J co u = Ir = 0. Only I a ^ 0. Furthermore 
S A (X,p) = S A (A) = m(X) alone, so that 

I A = Y,\X)[l- l ^(dJ p -*dJ x )\S- + 

- s- + [\-^Cd x ^ p -'d p 'n x )]i: A {x) + o{h 2 ) 

(141) 

or 

I A = -ihd^ A d£S- + , (142) 
and the transport equation becomes 

|{7 ^% + ^ P) }+[^^ + (A, P )] = -iW^ p S-\ (143) 



Assuming that the quasiparticle ansatz for S + (X,p) of Eq.( |140|) holds and 



that the mass is to be considered as the dynamically generated Hartree 
mass that is to be self-consistently determined, one can insert Eq. (|140|) into 
Eq.( |143j ), take the trace over spinor indices and integrate over a positive en- 
ergy interval A + that contains E p , to arrive at an equation for the quark 
distribution function, 

^f q (X,p) = Pid J -^- + m(X)dMX)d;(^^) = 0. (144) 
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On performing the derivatives and extracting a factor of l/E p , one can write 
this as 

-L(^/ ? (X,£) + m(X)d fl m( X Hf q {X,p)) = 0, (145) 
tip 

which is the Vlasov equation for the model. It must be solved concurrently 
with the gap equation for m(X), 

m(X) = m + AGN c m(X) J * [1 - f q (X,p) - f q (X,p)}, (146) 

that is derived directly from the Hartree self-energy. The constraint equation, 
in this same approximation in the expansion in h, is 

{p 2 -m 2 {X))f q (X,p) = 0, (147) 

which validates our use of the quasiparticle assumption as being exact. Equa- 
tion (|145|) indiates that chiral symmetry breaking enters via the condensate 



or mass already as a spatially varying potential in the Vlasov equation. 

4.4 The Boltzmann equation for the NJL model 

In principle, the next step from a physical point of view would be to in- 
corporate all self-energy diagrams of the next order in 1/N C . This would 
correspond to meson exchange ||12|| . This has not been done yet formally 
|4"?f and we will touch on this briefly in the following subsection. Here we 
shall rather examine the simpler problem of considering our self-energy with 
at least two interaction vertices, such as shown in Figs. 17 and 18 for the 
NJL model. These are the minimal types of diagram that can possibly give 
rise to an off-diagonal self-energy £ + ~ say, and therefore to a non-vanishing 
contribution to the gain and loss terms that comprise 7 co n in the transport 
equation, Eq. (|133|) . We will not give details here, but just note the salient 
features |J46|. 



52 



(a) (b) 

Figure 17: Direct (a) and exchange (b) graphs that contribute to E~ + and 
which contain two interaction lines. The vertices can be either all scalar or 
all pseudoscalar in nature. 




Figure 18: Mixed graphs that contribute to E^ ixcd and which contain two 
interaction lines. 



Firstly, a direct translation of the off-diagonal graphs, in the scalar chan- 
nel say, 

K (X,P) = -4G h J ^ ( ^ )4 (2irft) 5(p - Pl+P2 - p 3 ) 

x [S + -(X, Pl )tr (S- + (X,p 2 )S + -(X,p 3 )) 

(X, Pl ) S~ + (X, p 2 ) S + ~ (X,p 3 )}, ( 148) 

contains a product of three Green functions. Recalling that this will be 
multiplied by S~ + in the collision integral and also that we must trace and 
integrate the result first over a positive energy interval, we can easily see that 
such a procedure will lead to eight terms that each contain some product of 
four quark or antiquark distribution functions, such as for example 

coefficient x f q {p)f q {p2)f q {pz)f q {pi) (149) 

In a loose sense, if one designates f q ,q(p) to represent an incoming quark 
(antiquark) and f qq to represent an outgoing quark (antiquark), then one 
can draw diagrams associated with each process. For example, the products 
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(a) (b) 




(c) (d) 

Figure 19: Six graphs that arize from the term ~ H + ~S h . These are 
heuristic graphs are are not Feynman diagrams. 

listed in Eq. ( |149| ) would represent quark-quark scattering. A similar term of 
the eight possible leads to quark- ant iquark scattering, while the remaining 
six that are not listed (but which are easily worked out), are shown in Fig. 19. 
Some of these look like the typical vacuum fluctuation processes that would 
occur in any relativistic theory and in addition to these, there are others 
that give rise to pair creation and annihilation. All six graphs of this figure 
can be shown to vanish from energy-momentum conservation due to the 
quasiparticle assumption! This gives us an indication of the complexity and 
richness of the theory that would go beyond the standard collision scenario 
if one relaxes this assumption. 



Secondly, it is important to verify that the coefficient functions in the 
term of (|149|) in fact truly give rise to the differential cross section for elastic 
quark-quark scattering as would be calculated from real Feynman diagrams 
(and not heuristic graphs of Fig. 16) such as are displayed in Fig. 20. In 
fact, this has been explicitly demonstrated to be the case [jlql . One finds 
that the contribution from Fig. 17(a) gives rise to the amplitudes squared of 
both the s or u channels for qq scattering (or s or i channels for qq scatter- 
ing), while Fig. 17(b) is required to produce the interference terms between 
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P 2 P 2 

^-channel li-channel 



Figure 20: t and u channel Feynman graphs for elastic quark-quark scatter- 
ing, to lowest order in 1/N C . 



them. It appears that evaluating nonequilibrium self-energies for the Boltz- 
mann equation leads to scatering processes that can be obtained from all 
possible combinations of cutting the slef-energy grphas of Fig. 17 vertically, 
reminiscent of the Wick-Cutkowsky rules ||17||. 



Finally, one arrives at a Boltzmann equation from Eq. (|133|) . It reads 



P"SMX,P) + m(X)d„m(X)d£f,(X,p) 



l>2 



l>2 



+ ^I M --W-(P2 -> 13)(/ g (pi)/ ? (p2)/ff(P3)/ 9 (p) - fq(Pl)fq(P3)fq(P2)fq(p))}, 



(150) 



The constraint derived earlier, Eq. (|147|) , however, remains unaltered. From 
the Boltzmann equation, it is apparent that the changes in the condensate 
with the medium affect the equation in two possible places: (a) As with 
the Vlasov equation, a medium dependent potential occurs on the left hand 
side that is related to the effective quark mass in medium and (b) the cross- 
sections occurring on the right hand side are medium dependent, and also 
depend on changes of the quark and meson masses in the medium. As we 
have seen in the preceding section, the cross section for quark-antiquark 
scattering diverges at the phase transition. 
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Figure 21: Self-consistent self-energy that includes meson exchange. 



The actual answer as to what one should expect from numerical simu- 
lations is however unclear: since the differential cross-sections are averaged 
over, one may lose the sharp signal of the divergence. Howver, the force term 
on the left hand side may still play an essential role. At this stage also, too 
many physical features are still lacking, in particular, the coupling of the 
quark degrees of freedom to mesons and their coupling back to the quarks. 
This must lead to a hadronization scenario. In the final subsection of this 
chapter, we briefly sketch how this might occur. For numerical simulations 
thus far, we refer the reader to [48] and other references cited therein. 



4.5 Higher orders in 1/N C and meson production 

As already pointed out earlier, the expansion in the coupling strength that 
was used for selecting the diagrams of the last section is inadmissable, because 
GA ~ 2. Going to higher orders in the 1/N C expansion is however non-trivial, 
as a symmetry conserving set of graphs must be chosen. From [O, [13|, we 



know that this comprises firstly the set of graphs of Fig. 21 for the self- 
energy, where the "F" denotes the new full Green function that must be 
newly determined in a self consistent fashion. 



Denoting the two terms in the self-energy as Si and E 2 respectively, one 
can make an expansion of the new full Green function about that governed 
by £1, 

S = S 1 + S 2 (A;) (151) 
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Figure 22: Contributions to the irreducible polarization to next to leading 
order in the 1/N C expansion. 



= + < 152 > 

Concomitantly, the irreducible polarization U F (k) now occurring in the quark 
antiquark scattering amplitude 

- wF ^-T^m (153) 

must contain further terms, 

n F = n° + <m (154) 

where II is the simple quark loop, in order to be symmetry conserving. The 
graphs required for U F are shown in Fig. 22. 



Inserting the expansion of the Green function and the irreducible polar- 
ization into the full self-energy of Fig. 21, leads to graphs that contain inter 
alia diagrams of the form shown in Fig. 23. 
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Figure 23: Diagrams occurring in E, that lead to the hadronization of a 
quark into two mesons. 



This gives us an intuitive understanding that, on evaluating these dia- 
grams in the non-equilibrium scenario, we should no longer simply obtain 
a cross-section for elastic quark-quark and quark- ant iquark scattering, but 
also the hadronization process of qq — > MM', where M and M' are mesons. 
Much work however, remains to be done in this regard. 

5 Concluding comments 

In this series of lectures, we have investigated some aspects of chiral symmetry 
breaking at finite temperatures. We have seen that in the last few years, 
much information is emerging from the lattice gauge community that tells us 
about the transition region itself. Chiral perturbation theory, on the other 
hand, while being excellent in the low temperature regime, cannot adequately 
describe a phase transition. 

In the following section, we have investigated the Nambu-Jona-Lasinio 
model at finite temperatures. It gives a remarkably good qualitative agree- 
ment with the lattice data in the realm of static properties. It fails, however, 
to describe the bulk thermodynamic properties well, primarily due to the 
fact that confinement is lacking. The NJL model gives a simple picture for a 
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derealization rather than a deconfinement transition. Associated with this 
(physically appealing) picture that bound mesons become delocalized at the 
transition temperature - now the Mott temperature - and are still correlated 
states with a finite width in the quark medium, are marked divergences in 
many functions, such as the pion radius, tt-tt and tt-K scattering lengths (not 
discussed here), as well as the phenomenon of critical scattering, observed in 
the quark- ant iquark channel. 

Due to the fact that none of the apparent singularities are directly ob- 
servable experimentally, we have turned to transport theory, in order to in- 
vestigate what effects are to be expected from a condensate density that is 
medium dependent. Calculations at this stage indicate that a Boltzmann 
equation is dependent on the condensate through a force term, and also via 
the cross-sections that arize from binary collisions among the quarks and 
antiquarks. Howver, the stage of calculation is still primitive: a consistent 
physical theory that includes mesons and which overcomes the problems as- 
sociated with the lack of confinement is required before one can expect to 
obtain credible results. This, of course, leaves the path open for future re- 
search. 
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